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A  new  remote  sensing  technique  was  developed  for  pre- 
dicting the  expected  mean  laser  radar  return  from  a  rough 
water  surface.  This  technique  involved  measuring  the  stan- 
dard deviations  of  the  upwind  and  crosswind  profiles  of  the 
elliptical  glitter  patterns  occurring  for  illumination  of 
the  water  surface  with  a  point  source  near  the  laser  radar 
system.  A  pencil  beam  laser  radar  from  a  companion  project 
simultaneously  measured  the  reflected  signals  from  the  water 
surface.  The  glitter  pattern  images  were  recorded  with  a 
video  camera  and  recorder.  The  images  for  each  run  were 
later  digitized  along  their  major  and  minor  elliptic  axes 
and  averaged  over  256  images  to  produce  smooth  intensity 
curves  from  which  the  standard  deviations  were  measured. 
The  radar  return  fluctuated  over  a  large  range  because  of 
the  rapid  variation  of  individual  water  surface  facets,  and 
so  was  recorded  and  time-averaged  over  the  same  interval  as 
the  video  images.  Data  sufficient  for  empirical  prediction 
of  expected  mean  laser  return  signal  were  obtained.  This  is 
necessary  to  permit  evaluation  of  the  performance  of  a  given 
laser  radar  design.  The  data  obtained  also  approximated  the 
predictions  of  a  new  model  proposed  in  this  work. 
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I.     INTRODUCTION 

A .      BACKGROUND 

The  situation  at  the  outset  of  this  project  was  that  a 
working  model  of  an  inexpensive,  light,  and  compact  laser 
radar  altimeter,  for  use  over  the  ocean  surface,  had  been 
developed  in  this  laboratory  in  a  companion  project 
[Ref.  1],  However,  it  had  not  been  possible  to  relate  its 
performance  in  the  laboratory  to  real  field  performance  with 
much  certainty.  The  device  determined  the  height  above  the 
water  by  measuring  the  time  delay  between  the  emission  of  a 
short  laser  pulse  and  the  arrival  of  a  return  echo  from  the 
water  surface.  The  limitations  of  size,  weight,  electrical 
power,  and  cost  dictated  a  low  power  diode  laser  system. 
Additionally,  there  was  a  requirement  of  a  finite  divergence 
angle  of  the  laser  beam  and  associated  detector  optics  in 
order  to  assure  triggering  of  the  time-delay  measuring  sys- 
tem in  the  event  of  a  glassy  calm  and  some  angular  wobble  in 
the  attitude  of  the  laser  system.  These  collective  require- 
ments meant  that  the  primary  design  problem  was  that  of 
achieving  sufficient  return  signal  from  the  water  to  assure 
triggering    of    the    time    delay    circuits. 

Evaluating  the  performance  of  the  existing  working  model 
of  a  laser  radar  altimeter  system  over  the  water  had  the 
problem  that  the  reflectance  of  a  water  surface  is  dependent 
on  the  roughness  state  of  the  water.  The  only  data  in  the 
literature     known     to     this     author,      which     would     permit     an 


estimate  of  the  laser  reflectance  of  rough  water,  are  those 
of  Petri  [Ref.  2],  These  data  are  rather  meager,  consis- 
ting of  16  separate  measurements,  and  requiring  a  know- 
ledge of  the  wind  speed.  The  latter  had  been  measured  at  a 
height  of  60  feet  above  the  water.  The  wind  speed  at  the 
water  surface  was  thus  quite  uncertain.  Other  data  by  Cox 
and  Munk  [Ref.  3]  related  the  glitter  pattern  profile,  for 
reflection  of  the  sun,  to  wind  speed  over  the  ocean.  Their 
optical  data  were  taken  from  an  aircraft  at  2000  feet  alti- 
tude, with  windspeed  measured  on  a  ship  at  two  heights,  9 
feet  and  41  feet  above  the  water.  It  was  not  clear  which 
of  these  two  heights  was  used  for  quoting  the  wind  speed. 
Their  data  also  were  not  directed  at  evaluating  the  magni- 
tude of  the  reflectance. 

B.   GENERAL   OBJECTIVE 

Preliminary  field  experiments,  carried  out  as  part  of 
the  work  reported  here,  indicated  that  the  reflectance  of 
rough  water  could  vary  rapidly  under  changing  wind  condi- 
tions, and  that  knowledge  of  the  wind  speed  at  a  given 
instant  did  not  serve  as  a  good  indicator  of  the  optical 
properties  of  the  water  surface  at  that  instant.  Addi- 
tionally, if  wind  speed  were  to  be  measured,  it  probably 
should  be  measured  very  close  to  the  water  surface.  This 
would  be  difficult,  expensive,  and  in  many  cases  impracti- 
cal. The  field  work  reported  here  for  this  project  was 
carried  out  from  high  bridges,  but  it  was  intended  that  the 


work  would  be  extended  later  to  aircraft.  In  that  case  the 
necessity  for  associated  ship  measurements  of  wind  speed 
would  be  difficult  logistically.  Consequently,  the  objec- 
tive of  this  work  became  to  develop  and  demonstrate  the 
feasibility  of  a  remote  sensing-  technique  for  evaluating 
the  reflectance  of  a  rough  water  surface  at  the  instant  of 
a  laser  radar  altimeter  test. 

The  papers  of  Cox  and  Munk,  and  Petri,  taken  together, 
led  to  a  conclusion  that  it  might  be  possible  to  predict  an 
expected  laser  radar  return  on  the  basis  of  a  measurement 
of  the  width  of  glitter  patterns.  It  also  seemed  likely 
that  glitter  patterns  measured  with  a  nearly  point  source 
of  light  at  night  would  be  simpler  to  obtain  and  would 
yield  more  directly  the  desired  result  than  would  glitter 
patterns  obtained  with  the  necessarily  oblique  illumination 
from  the  sun.  Consequently,  the  primary  objective  became 
to  test  the  feasibility  of  prediction  of  laser  return  from 
the  profile  widths  of  glitter  patterns  obtained  with  a 
point  light  source.  The  early  tests  indicated  that  laser 
return  was  in  fact  a  function  of  glitter  pattern  width.  A 
tentative  theoretical  model  was  developed,  which,  although 
it  has  not  as  yet  succeeded  in  predicting  the  observed 
absolute  magnitude  of  laser  return,  does  lead  to  a  func- 
tional relationship  which  is  useful.  With  the  present 
semiempir ical  relationship,  laser  return  can  be  related 
quantitatively  to  Lambertian  return.  Although  more  experi- 
mental data  is  needed  to  refine  the  results,  and  further 
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theoretical  work  is  needed  to  clear  up  the  discrepancy  in 
the  absolute  values  predicted  by  the  model,  laser  return 
signals  measured  in  the  laboratory  with  Lambertian  targets 
can  now  be  directly  related  to  field-expected  values. 

C.   ASSOCIATED  PROBLEMS 

1.   Laser  Signal  Fluctuation 

There  were  several  additional  problems  in  carrying 
out  the  field  experiments  at  the  outset  of  this  project. 
The  laser  radar  altimeter  model  to  be  tested  used  a  fairly 
narrow  laser  beam  of  1.2  by  2.4  degrees  divergence,  with  a 
3  degree  circular  cone  of  acceptance  of  the  receiver  op- 
tics. This  divergence  had  been  chosen  in  the  design  pro- 
cess of  the  companion  project  as  an  optimization  of  the 
transmitter-receiver  optics  under  the  confines  of  required 
return  signal  magnitude  and  availability  of  inexpensive 
commercial  components.  This  narrow  beam  was  incident  on 
the  center  of  the  "glitter"  pattern,  about  which  much  more 
will  be  said  later.  The  beam  will  be  treated  here  as  a 
pencil  beam,  essentially  evaluating  the  peak  value  of  laser 
return  at  the  center  of  the  glitter  pattern.  A  consequence 
of  the  narrow  beam  geometry  was  that  reflection  of  the 
laser  beam  from  the  water  surface  occurred  at  a  relatively 
small  and  finite  number  of  glint  facets.  As  a  result,  the 
laser  return  signal  fluctuated  rapidly  with  time  over  a 
wide  range  of  magnitude.  This  required  recording  the 
signal  for  later  data  processing  to  obtain  the  mean  and 
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standard  deviation  of  the  variation.  In  order  to  accom- 
plish the  recording  it  was  also  necessary  to  provide  a 
pulse  stretcher  because  the  laser  return  pulses  were  other- 
wise too  short  to  handle  with  most  analog  FM  data  recorders 

Although  the  rapidly  fluctuating  magnitude  of  the 
laser  return  posed  a  measurement  problem,  it  is  helpful  in 
the  practical  application  of  a  laser  radar  altimeter.  Mea- 
surement of  range  is  accomplished  by  determining  the  time 
delay  between  the  primary  pulse  and  the  return  pulse.  With 
fairly  rapid  pulse  repetition  rate,  an  occasional  high 
pulse  will  serve  to  make  the  range  determination  practical, 
where  the  mean  signal  may  not  be  above  threshold  for  satis- 
factory operation.  This  feature  tends  to  favor  a  small 
divergence  beam.  The  statistics  of  this  were  not  analyzed 
here  but  are  left  for  later  work.  The  data  tape  recordings 
from  the  field  experiments  are  now  available  for  such 
analys  is  . 

2 .   Calibration  by  Means  of  Specular  Reflection 

Prior  to  the  work  reported  here,  optimization  of 
design  and  laboratory  evaluation  of  the  laser  radar  altime- 
ter had  been  carried  out  primarily  using  specular  reflec- 
tion from  the  outside  of  a  plate  glass  window  to  obtain  a 
return  signal,  after  careful  measurement  of  the  inevitable 
slight  surface  curvature  of  the  window.  This  permitted 
development  of  the  time-difference  circuitry  and  avoided 
the  difficult  problem  of  obtaining  a  water  surface  at 
sufficient  range  to  be  meaningful.   The  constant  magnitude 
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of  the  return  signals  in  this  situation  was  also  helpful  in 
the  development  work.  In  evaluating  the  magnitude  of  the 
signal,  the  difference  of  reflectance  of  the  two  glass-air 
interfaces  from  that  of  a  single  water  surface  was  easily 
accounted  for.  A  rough  evaluation  of  system  performance 
was  obtained  in  this  manner,  but  the  return  signal,  if 
maximized  by  angular  adjustment  to  obtain  the  maximum  sig- 
nal, evaluated  only  the  peak  value  of  the  central  maximum 
of  the  laser  beam,  as  will  be  discussed  in  Appendix  B.  To 
obtain  a  measure  of  the  expected  return  signal  from  a  rough 
water  surface  or  a  diffuse  reflecting  surface  from  the 
specular  reflection  data  required  carrying  out  a  two-dimen- 
sional integration  of  the  laser  beam  flux  over  the  laser 
beam  cross-section  profile.  The  beam  tended  to  have  consi- 
derable detailed  "banding"  structure  in  its  profile,  making 
such    evaluation    rather    doubtful. 

3.        Calibration    by    Means    of    Lambertian    Reflection 

The  problems  attendant  on  use  of  specular  reflec- 
tion from  a  plate  glass  window  to  obtain  a  return  signal 
are  alleviated  by  use  of  a  diffuse,  "Lambertian"  reflector 
such  as  a  white  sheet  of  paper.  For  a  finite  area,  such  a 
surface  reflects  a  power  proportional  to  the  cosine  of  the 
angle  of  incidence.  Such  surfaces  were  used  in  the  labora- 
tory. However,  the  Lambertian  surface  must  be  large  enough 
to  include  the  entire  transmitted  laser  spot.  Hence  at 
ranges  sufficient  to  simulate  real  field  conditions,  the 
area    required    was    large    enough    to    make    obtaining    a    good 
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uniform  reflectance  white  surface  difficult.  The  signals 
with  such  a  target  would  be  constant  in  magnitude,  except 
for  a  small  amount  of  atmospheric  scintillation.  At  close 
range,  the  signals  were  so  large  as  to  make  circuit  devel- 
opment for  small  signal  conditions  difficult.  Calibration 
of  the  return  signal  magnitude  can  be  done  with  a  Lamber- 
tian  surface  at  close  range,  by  careful  attention  to  quan- 
titative attenuation  of  the  large  return  signal  that  would 
otherwise  saturate  the  receiver  circuits.  This  sort  of 
calibration  was  done  in  the  field  tests  to  be  reported 
here.  Finally  there  seemed  to  be  no  fully  satisfactory  way 
to  simulate  all  the  circumstances  to  be  encountered  by  a 
real  life  laser  radar  altimeter,  other  than  by  field  tests 
over  wind-roughened  water.  These  were  then  initiated  and 
the  accompanying  remote  sensing  techniques  reported  here 
carried  out. 

D.   GENERAL  METHOD 

This  work  reports  a  new  remote  sensing  system  for  esti- 
mating the  expected  laser  radar  return  signal  through  an- 
alysis of  glitter  patterns.  This  is  a  remote  sensing 
method  permitting  determination  of  the  expected  reflectance 
of  a  rough  water  surface  using  equipment  colocated  with  the 
laser  radar  altimeter.  The  equipment  for  this  technique 
also  offers  the  advantage  of  being  inexpensive  and  phy- 
sically small  for  field  use. 
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The  time-averaged  statistical  distribution  of  light 
intensity  was  measured  for  images  of  glitter  patterns  pro- 
duced by  reflection  of  light  projected  onto  the  water  from 
a  small  area  incandescent  source  located  near  the  laser 
altimeter.  The  images  of  the  glitter  patterns  were  re- 
corded with  a  video  camera,  which  avoided  the  nonlineari- 
ties  of  photographic  recording  and  permitted  digitization 
without  intermediate  steps.  The  magnitude  of  the  returned 
laser  signal  was  recorded  at  the  same  time  that  the  video 
images  were  obtained.  Because  this  signal  was  rapidly 
varying,  it  was  recorded  for  statistical  processing  later. 
The  resulting  data  related  the  mean  laser  return  to  the 
glitter  pattern  profile  parameters.  The  results  are  most 
succintly  expressed  as  a  ratio  of  reflectance  of  the  rough 
water  to  the  reflectance  of  a  (virtual)  Lambertian  reflec- 
tor which  reflects  the  same  fraction  of  the  total  light  as 
a  water  surface  (0.0204)  at  normal  incidence.  Through  use 
of  data  in  other  references,  the  return  was  also  related  to 
an  inferred  wind  velocity  for  comparison.  An  analytical 
model  was  developed  which  seems  to  predict  the  correct 
functional  behavior,  although  the  predicted  absolute  magni- 
tude of  reflectance  is  at  present  out  of  line.  Although 
these  experiments  should  be  repeated  for  refinement  and 
completion,  these  results  allow  evaluation  of  a  practical 
laser  radar  system  by  reflection  from  a  Lambertian  reflec- 
tor within  the  laboratory,  providing  a  more  immediate  and 
readily  controlled  environment  for  equipment  development. 
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II    THEORETICAL  CONSIDERATIONS 

A.    GENERAL 

When  light  from  a  small  diverging  source  is  incident  on 
a  water  surface  roughened  by  wind  driven  waves,  it  is  re- 
flected from  the  water  surface  in  a  glitter  pattern,  as 
seen  by  an  imager  located  near  the  source.  The  glitter 
pattern  is  approximately  elliptical  in  shape,  with  the 
major  axis  of  the  ellipse  in  the  upwind  direction,  and  the 
minor  axis  in  the  crosswind  direction  [Ref.  4],  The  pat- 
tern is  produced  by  specular  reflection  at  a  large  number 
of  rapidly  changing  facets  located  at  the  points  where  the 
surface  is  perpendicular  to  the  incident  light.  A  repre- 
sentative  glitter    pattern    is    shown    in    Figure   1. 

Laser  reflectance  from  calm  waters  can  be  readily  cal- 
culated analytically  [Ref.  5],  This  is  the  case  of  specu- 
lar reflection  whereby  a  mirror  image  of  the  source  is 
created  equidistant  below  the  water  surface.  Reflected 
energy   appears    to  originate   from   this    source. 

At  the  other  extreme  from  specular  reflection  is  dif- 
fuse reflection.  In  this  case  a  surface  scatters  incident 
energy  over  all  angles  in  a  hemisphere,  resulting  in  a 
distribution  called  Lambertian,  in  which  the  reflected 
energy  varies  as  the  cosine  of  the  angle  from  the  perpen- 
dicular. 


"< 
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Figure  1.   Glitter  pattern  produced  by  wind  speed 
of  approximately  12  miles  per  hour. 
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Between  these  two  extremes,  an  optically  rough  water 
surface  exhibits  properties  of  both  types  of  reflection. 
The  reflected  intensity  distribution  varies  according  to 
the  distribution  of  wave  slopes  scattering  the  incident 
beam.    [ Ref .  6  ]  . 
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B.      GLITTER   PATTERNS 

The  technique  developed  in  this  work  was  a  method  of 
evaluating  the  expected  laser  radar  return  from  a  water 
surface  when  it  is  roughened  by  wind.  The  technique  is  a 
remote  sensing  technique  in  that  all  measurements  are  made 
from  the  general  location  of  the  laser  radar  altimeter. 
Thus  it  does  not  involve  knowing  the  windspeed  at  the  water 
surface,  nor  any  other  meteorologic  or  oceanographic  par- 
ameters. The  measurements  involve  only  the  optical  proper- 
ties of  the  water  surface.  This  is  done  by  determining  the 
size  of  the  glitter  pattern  from  a  video  camera  recording 
of  the  reflection  of  an  approximately  point  source  of  light 
near  the  laser  radar.  The  laser  radar  used  a  pencil  beam, 
so  that  it  responded  to  the  central  maximum  reflection 
point  of  the  glitter  pattern  profile.  The  reflection  from 
this  fluctuated  through  a  large  range,  as  a  function  of 
time,  so  that  a  time-average  of  a  recorded  signal  was 
requi  red . 

Cox  and  Munk  have  developed  a  method  of  measuring  water 
surface  roughness  from  photographs  of  sun  glitter  [Ref.  3]. 
The  wave  slope  distribution  was  deduced  from  observations 
of  sun  glitter  patterns,  with  the  assumption  that  the  wave 
slope  distribution  and  the  reflected  light  intensity  dis- 
tribution were  the  same.  They  found  that  a  Gram-Char  1 ier 
distribution  was  the  best  fit  to  the  experimental  data. 
This    distribution    function    is 


19 


[l/(2nanO  ]    x   e 


-(X2+Y2)/2 


u   c 
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<-Skew 


Peak- 
edness 


+    (1/4)    C99(X2-1)(Y2-1)    +    ...] 


or,      G   =       [l/(2nau<rc)]    x    e 


x  [H] 


(1) 


where 


X  =  9/ffu   with  0  the  angle  from  vertical  in  the 
upwind/downwind  direction,  and 

Y  =  G>/<*c   ,  with  0  the  angle  from  the  vertical 
in  the  crosswind  direction. 

The  quantities  <*u  and  oQ  are  the  upwind  and  crosswind 
"standard  deviations"  for  this  distribution.  For  a  Gaus- 
sian distribution  the  standard  deviation  is  the  root  mean 
square  deviation  from  the  mean.  The  first  two  terms  of  the 
Gram-Char  1 i er  distribution  are  the  same  as  for  a  Gaussian 
distribution.  The  Gram-Char  1 ier  sigma  values  are  the  val- 
ues of  sigma  in  the  second  term.  As  will  be  discussed 
later,  the  peakedness  of  the  Gram-Char  1 i er  distribution 
causes  these  sigma  values  to  be  located  lower  on  the  curve 
than  for  a  Gaussian  distribution.  The  terms  in  the  brac- 
kets in  the  Gram-Char  1  ier  distribution  are  a  form  of  Her- 
mite  polynomial.  These  terms  express  the  peakedness  and 
skew  of  the  function. 

The  Gram-Char  1 i er  distribution,  above,  is  a  two  di- 
mensional distribution.   It  is  displayed  as  a  function  of 
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0  and  0  in  Figure  2.  The  dotted  curves  are  the  shape  of  a 
Gaussian  with  the  same  standard  deviation  values.  The 
crosswind  distribution  is  symmetric  about  the  vertical  but 
is  more  peaked  than  Gaussian.  The  upwind  distribution  is 
also  more  peaked  than  Gaussian,  and  has  a  skew  in  the 
upwind  direction. 

The  Gram-Char  1 i er  distribution  is  a  normalized  func- 
tion, i.e.,  with  the  complete  Hermite  polynomial  series 
included,  the  integral  of  the  function  from  minus  infinity 
to  plus  infinity  in  X,  and  from  minus  infinity  to  plus 
infinity  in  Y,  gives  unity.  At  first  sight  it  would  seem 
that  integration  from  minus  infinity  to  plus  infinity,  over 
the  angular  variables  used,  is  not  strictly  correct,  as  the 
angles  can  only  go  to  n/2.  The  function  as  written  above 
is  an  approximation,  as  it  has  only  the  first  few  terms  of 
the  Hermite  polynomial  included.  However,  Cox  and  Munk 
stated  that  this  approximation  was  valid  for  all  their  data 
out  to  X  =  Y  =  2.5.  Contributions  to  the  integral  for 
values  of  X  and  Y  beyond  2.5  are  negligible,  so  that  writ- 
ing the  integral  as  if  it  extended  to  infinity  should  be  of 
no  concern.  The  largest  sigma  values  encountered  in  the 
data  reported  here  were  7.5  degrees.  For  X  =  2.5  this 
means  the  expression  is  taken  to  be  valid  to  19  degrees. 
The  integrals  written  later  in  this  thesis  will  express  the 
limits  as  minus  infinity  to  plus  infinity,  recognizing  that 
the  contributions  of  high  X  and  Y  are  negligible. 
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Figure  2.  Wave  slope  distribution  profiles  deduced  by 
Cox  and  Munk  [Ref.  3],  Upper  curves  are  taken  along  the 
Minor  (Crosswind)  axis  of  the  elliptical  pattern.  Lower 
curves  are  taken  along  the  major  (up/downwind)  axis.  Solid 
curves  are  distributions  deduced  for  a  wind  speed  of  10 
meters  per  second.  Dashed  curves  refer  to  Gaussian  distri- 
butions   of    the   same   standard   deviation. 
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The  extent  to  which  the  Gram-Char  1  i  er  distribution 
exceeds  the  Gaussian  due  to  the  peakedness  at  the  center 
can  be  evaluated  by  inserting  the  values  of  the  coef- 
ficients, as  given  by  Cox  and  Munk,  and  letting  X  =  0  and  Y 
=  0.  For  the  crosswind  direction,  the  curve  has  no  skew 
and  the  peakedness  is  independent  of  wind  velocity.  The 
coefficients  for  crosswind  peakedness  are 
c40    =    0.40    ±    0.23    , 


04 


=    0.23    ±    0.41     , 


(2) 


and  c22    =    0.12    ±    0.06    . 

This  gives  a  term  multiplying  the  value  of  the  peak  as 
given  by  the  first  two  terms  of  equation  (1),  i.e.  the 
height  of  the  Gaussian  central  peak  is  multiplied  by  a 
factor  f    =    1.109    ±    0.061     .  (3) 

This    factor   will    be   used   as    1.11    later. 

For    the    upwind   direction    the    skewness    constants    are 
c21    =    0.01    -    0.0086W   ±    0.03    , 
and  c03    =    0.04    -      0.033W   ±    0.12    ,  (4) 

where  W  is  the  wind  speed  in  meters  per  second.  Although 
these  constants  depend  on  wind  speed,  they  do  not  affect 
the  area  under  the  curve,  and  hence  the  normalization,  as  X 
and  Y  appear  only  in  odd  powers  in  the  Gram-Char  1  ier  Hermi- 
tian  terms.  Thus  whatever  area  they  add  for  positive 
values  of    X   or    Y,    they   subtract    for    negative   X   or    Y. 

Although  the  skewness  terms  do  not  affect  the  peaked- 
ness for  the  crosswind  direction,  they  do  make  a  small 
contribution     to     the     peakedness     in     the     up/downwind 
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direction,  but  the  contribution  is  small.  For  example,  for 
a  wind  speed  of  10  meters  per  second,  as  illustrated  in  the 
lower  curve  of  Figure  2,  the  additional  peakedness  is 
approximately  two  percent.  The  position  of  the  peak  is 
also  shifted  slightly  upwind.  However,  its  position  is 
still  within  the  "pencil"  beam  of  1.4  x  2.8  degrees  used  in 
the  laser  radar  altimeter  here.  The  values  for  the  predic- 
ted laser  return  signal  presented  later  in  this  thesis  have 
been  adjusted  for  the  skewness  contribution  to  the  peaked- 
ness, even  though  this  adjustment  is  in  most  cases  the 
order  of  one  percent  -  really  negligible  in  view  of  the  6 
percent  uncertainty  quoted  above  from  Cox  and  Munk  on  the 
value   of    the   peakedness    factor. 

To  return  to  the  meaning  of  the  "standard  deviations" 
for  the  Gram-Char  1 ier  distribution,  the  standard  deviation 
values,  oQ  and  <xu,  are  the  standard  deviations  of  the 
Gaussian  term  of  the  Gr am -Cha r  1  i e r  form.  In  measuring 
these  quantities  on  an  actual  measured  distribution,  as 
will  be  described  later  in  this  work,  it  is  necessary  to 
measure  the  distribution  curve  width  at  a  point  lower  by  a 
fraction  1.11  (or  1.12,  depending  on  the  skew  contribution 
to  the  peakedness  for  that  particular  curve)  than  the  usual 
0.659  times  the  peak  value  -  the  case  for  the  usual  Gauss- 
ian distribution.  This  has  been  done  in  evaluating  the 
data    later    in    this    thesis. 

A  further  fine  tuning  should  also  be  noted.  The  peak- 
edness   terms    in    the    Gram-Char  1  ier    distribution    do    produce 
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upward,  or  downward,  displacements  of  the  curve  throughout 
the  curve.  This  can  be  seen  in  the  curves  presented  in 
Figure  2.  Because  of  this  displacement,  it  might  seem  that 
determination  of  sigma  values  by  measuring  the  width  of  the 
peak  at  a  height  equal  to  0.659/1.11  times  the  peak  value, 
would  be  in  error.  (The  factor  1.11  comes  from  equation 
(3)  above.)  To  evaluate  this  effect,  the  coefficients  in 
equation  (2)  above  produce  vertical  displacements  of  the 
crosswind  curve,  i.e.  for  Y  =  0  ,  at  X  =  +1  (at  the  o~c 
point) ,   given  by 

1    +    (1/24)    c40    (1-6+3) 
+    (1/24)    c04    (3) 

=   1   -    .0333   ±    .0192 
+    .0288   ±    .0513 

=  1  -  .0045  ±  .0547  . 
This  represents  a  downward  shift  of  0.4  percent,  with  an 
uncertainty  of  5.5  percent.  The  0.4  percent  correction 
could  be  made  by  drawing  a  parallel  line  to  the  curve, 
shifted  downward  by  that  fraction,  but  it  seems  not  worth 
doing  in  view  of  the  5.5  percent  uncertainty.  For  the 
up/downwind  curve  the  corrections  are  equally  insignificant 
Consequently  such  corrections  have  not  been  made  during 
data   reduction. 

The  most  significant  quantity,  for  this  work,  is  the 
ratio  of  the  peak  of  the  distribution  to  the  area  under  the 
distribution,  as  a  function  of  the  upwind  and  crosswind 
standard  deviations.      As  mentioned  before,    the  quantity,    G, 
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in  equation  (1)  above,  is  a  statistical  distribution.  It 
is  normalized  so  that  the  area  under  the  curve  is  equal  to 
unity,         i.e. 

Area    =    1    =J    J      [  1/ ( 2n<ruac)  ]    e  [H]    dX   dY ,       (5) 

where  the  Hermitian  peakedness  and  skewness  terms  have-  been 
abbreviated    by    the    bracket     [H], 

The  power  per  unit  solid  angle,  I,  reflected  from  a 
small  spot  on  the  water  surface  that  is  illuminated  by  the 
laser  radar  beam  will  have  its  reflected  light  spread  into 
a  cone  where  the  intensity  distribution  varies  according  to 
the  Gram-Char  1  ier  distribution.  The  magnitude  of  the  solid 
angle  of  the  cone  will  be  taken  to  be  that  subtended  by  the 
glitter  pattern,  as  observed  at  the  detector.  This  assump- 
tion   is   discussed   in  Appendix   D. 

Letting  IQ  be  the  power  per  unit  solid  angle  reflected 
perpendicular  to  the  water  surface  and  hence  back  to  the 
detector  of  the  laser  radar,  the  power  per  unit  solid  angle 
in   any   direction   will    be   given    by 

I    =    IQ    {2nauac/[H0]}    G    ,  (6) 

where    [HQ]    =    1.11    for    peakedness    at    X  =    0    and  Y  =    0    . 

As  a  verification  of  the  statement  in  equation  (6), 
note  that  at   X   =  0   and   Y   =   0    , 

G  =   GQ    =    {l/(2n<xuac)}    [HQ  ]    . 
Hence   at    X  =    0    and  Y  =    0    , 

I    =    I0    {(2nau<rc)/[H0]}     {l/(2nauac)}     [HQ]     , 
and         I    =    IQ    ,    as    it    should    be. 

Now  the    total    reflected   power    is    given   by 
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•  CO  /-CO 

p  = 


-oo     -co 


Jew    —  oo 
I    dX  dY    . 

Inserting    I    from  equation    (6)    gives 

lo    <2n<ru<xc/[H0])    G  dX  dY    . 

Then  from  equation  (5) 

P  =  I0  {(2Ttcruac/[H0]}  .  (7) 

Now,  the  total  reflected  power,  P,  is  the  power  radiated  by 

the  lasar  radar  in  the  pencil  beam,  PQ,  multiplied  by  the 

Fresnel  reflectance  of  water  at  normal  incidence,  r  ,  or 

P  =  rPQ  =  rl0  {2nauac/[H0]}  .  (8) 

Inverting  gives 

IQ    =    rPQ    (l/2nau<rc)     [Hj     .  (9) 

The    flux  at    the   detector,    FQff,    the    power    returned    per    unit 

area   of    the    detector    receiving   optics,    is    given    by    IQ    times 

the    solid    angle    subtended    by    a    unit    area    at     the    detector 

opt  ics ,    or 

Fog   =    I0   x    1/R2 

=  rPQ    {l/(2n0uac)  }    [HQ]   /R2   .  (10) 

C.      FRESNEL  REFLECTANCE 

Use  of  the  water  reflectance  for  normal  incidence  in 
the  preceding  section  perhaps  needs  some  discussion.  Dur- 
ing the  measurement  of  a  glitter  pattern  the  reflection  at 
each  glint  is  at  normal  incidence.  Similarly,  for  the 
light  scattered  back  to  the  laser  detector  from  a  pencil 
laser  radar  beam,  the  angle  of  incidence  is  again  zero.  On 
the  other  hand,  for  the  light  that  does  not  return  to  the 
receiver    because    it    is    reflected   out    into    an    elliptic    cone, 
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the  angle  of  incidence  at  the  reflecting  facets  is  not 
zero.  However,  the  contribution  of  scattered  light  for 
facets  with  angles  of  incidence  greater  than  2  sigma  is 
negligible.  Since  the  maximum  sigma  encountered  in  this 
work  was  7.5  degrees,  the  variation  of  reflectance  for 
angles  of  incidence  of  this  magnitude  needs  to  be  examined. 
The  reflectance  of  a  dielectric  surface  of  index  of 
refraction,  n,  as  given  by  the  Fresnel  equations,  for  light 
with  its  electric  vector  perpendicular  to  the  plane  of 
incidence    is 

cos  0  -   (n2  -  sin2  Q)1^2     2 


rs  = 


1/2 


(11) 


cos  0  +  (n2  -  sinz  0) 
For    light   with    its    electric   vector    parallel    to    the    plane   of 
incidence,    the    reflectance    is 

-n2cos  0  +  (n2  -  sin2  0)1/2 


rP  = 


1/2 


(12) 


rp  "  r» 


n^cos  0  +  (n2  -  sinz  0) 
At  normal  incidence,  0  =  0  and  r 
For  unpolarized  light  at  other  angles,  the  total  reflec- 
tance is  (rs  +  rp)/2  .  The  reflectance,  r,  at  some  repre- 
sentative angles,  for  water  of  index  n  =  1.333,  at  16C,  is 


e       rs      rD     r 

Normal  incidence     0°  .02037  .02037  .02037 

10°  .02133  .01944  .02038 

20°  .02449  .01662  .02055 

30°  .03093  .01194  .02144 

40°  .04316  .00585  .02450 

45°  .05299  .00281  .02790 

Brewster's  angle    53.1°  .07817  .00000  .07817 


Factor 

1.000 

1.001 

1.009 

1.053 

1.203 

1.370 

3.838 


Grazing  incidence   90 


1.0 


As  can  be  seen  from  the  above,  the  total  reflectance  for 
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unpolarized  light  varies  very  slowly  with  angle  of 
incidence,  even  up  to  fairly  large  angles.  The  increase  in 
reflectance  is  only  5%  at  30  degrees.  This  is  because  rp 
falls   as   rs   rises    for   angles    less    than    the    Brewster   angle. 

The  results  of  Cox  and  Munk  were  also  obtained  using 
the  same  assumption  as  made  here,  i.e.  that  the  reflectance 
of  the  water  surface  was  constant  and  equal  to  the  normal 
incidence   value    for   all    the    angles    encountered. 
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D.       LAMBERTIAN  REFLECTION 

In  order  to  compare  calculations  of  predicted  returned 
power  to  actual  measured  values,  it  is  useful  to  relate 
glitter  reflection  to  Lambertian  reflection,  where  scatter- 
ing is  completely  diffuse,  as  from  a  white  sheet  of  paper. 
However  the  comparison  will  be  made  to  a  Lambertian  type 
reflecting  surface  which  reflects  the  same  fraction  of  the 
total  power  on  it  as  the  reflectance  for  normal  incidence 
from  water,  i.e.  a  fraction  .0204  of  the  incident  light  is 
returned,  when  integrated  over  all  angles  within  a  hemis- 
phere . 

In  general,  for  Lambertian  reflection,  the  power  per 
unit  solid  angle,  I0l»  returned  at  perpendicular  incidence, 
is    given    by 

IoL   =   P/n  (13) 

where    P    is    the    total    power    reflected    from    the    surface, 
integrated   over   all    angles. 

In  the  case  of  an  idealized  Lambertian  water  surface, 
the   reflected   power   is   given   by 

P  =  rP0 
where,    as    before,    r    =    .0204    for    the   water    surface. 

The  flux  at  the  detector,  the  power  returned  per  unit 
area,    would  then  be  given  by 

FoL    =    rP0/"R2  •  <14> 

This   will   be  called   the   "Idealized   Lambertian    return    flux". 

To   summarize,    it    is    the   flux   returning  at    normal    incidence 
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from  a  Lambertian  surface  that  reflects,  integrated  over 
all  angles,  a  total  fraction  of  the  incident  light  equal  to 
the  reflectance  of  water  for  normal  incidence  (.0204). 

E.   EFFECTIVE  REFLECTANCE 

In  order  to  characterize  the  optical  state  of  the  rough 
water  surface,  it  is  useful  to  define  a  quantity,  Peff>  the 
effective  reflectance  of  the  water  surface,   as 

Return  flux  from  glitter 

Peff  =  

Idealized  Lambertian  return  flux 


F 
og 


(15) 


FoL 


This  quantity  can  be  measured  directly  in  terms  of  the 
ratio  of  laser  radar  return  signals  from  a  water  surface  to 
the  signals  from  a  Lambertian  surface  such  as  a  white  paper 
target,  often  at  a  different  range  from  that  of  the  water. 
This  is  discussed  in  detail  later  in  this  thesis,  but  to 
avoid  confusion  it  should  be  mentioned  here  that  the  ranges 
of  the  water  and  the  Lambertian  target  enter  in  that  case, 
as  well  as  the  actual  reflectance  of  the  Lambertian  surface 
used  as  a  reference.  In  turn,  the  anticipated  signal  from 
a  given  model  of  a  laser  radar  altimeter  can  be  estimated, 
using  a  knowledge  of  the  probable  range  of  Peff  to  be 
expected . 

In  terms  of  the  model  discussed  in  section  A  above,  an 
idealized,  or  predicted,  value  of  Peff  would  be  given  by 
equation    (15)    after    substituting    expressions    (10)    and    (14), 
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i  .e . 


or 


peff    = 


peff    = 


og 


oL 


[H0] 


2auac 


(16) 


This  will  be  called  the  idealized  Peff.  The  quantity  [  HQ  ] 
is  the  magnitude  of  the  Hermite  polynomial  for  the  two- 
dimensional  Gram-Char  1  ier  distribution  at  normal  incidence. 
[  HQ  ]  had  the  numerical  value  of  1.11  in  all  cases  encoun- 
tered here.  The  quantity  Peff  expresses  the  optical  pro- 
perties of  a  rough  water  surface,  in  terms  of  the  standard 
deviations,  <*u  and  <xc,  of  the  glitter  pattern.  It  is 
particularly  useful  because  it  does  not  depend  on  range  and 
depends  only  on  the  optical  state  of  the  rough  water  sur- 
face . 

The  functional  relationship  expressed  by  equation  (16) 
made  it  seem  desirable  to  plot  directly  measured  values  of 
Peff  as  a  function  of  l/(2<ru<y  c) .  The  results  of  such  a 
plot  are  discussed  later  in  this  thesis.  The  measured 
values    of    Peff    did    prove    to    vary    linearly   with    1/(2<tu<jc)  . 
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III.  EXPERIMENTAL  PROCEDURE 

A.   FIELD  SITE 

Experiments  to  measure  laser  reflectance  as  a  function 
of  glitter  intensity  profile  were  conducted  from  Parrott's 
Ferry  Bridge  over  the  New  Melones  reservoir  near  Columbia, 
California,  following  several  earlier  preliminary  experi- 
ments at  that  site  and  at  the  Dumbarton  bridge  over  San 
Francisco  Bay.  The  Parrott's  Ferry  Bridge  was  the  most 
favorable  of  many  sites  considered  as  it  offered  a  fairly 
large  distance  to  the  water  (109  feet  at  the  time  of  the 
last  experiments)  and  had  a  pedestrian  walkway  with  no 
obstructions  beneath.  Heavy  automobile  traffic  and  inade- 
quate pedestrian  space,  as  well  as  obstructions  beneath, 
made  such  sites  as  the  Golden  Gate  Bridge  and  other  bridges 
in  the  San  Francisco  Bay  area  unsuitable.  It  is  hoped  that 
funding  for  logistic  support  will  be  forthcoming  so  that 
these  experiments  can  be  continued  from  an  aerial  platform 
over  open  ocean. 

The  data  reported  here  were  all  collected  at  the  Par- 
rott's Ferry  Bridge  on  the  night  of  24  November,  1986,  a 
clear  night  with  no  moon  present.  The  wind  speed  varied 
from  approximately  3  miles  per  hour  to  12  miles  per  hour. 
Wind  direction  varied  less  than  30  degrees.  Due  to  the 
short  fetch  in  the  reservoir,  varying  from  1/4  mile  to  1 
mile,  depending  on  the  direction,  only  high  frequency  waves 
and    ripples    were    observed.       It    is    anticipated    that    results 
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with    these   waves    will    scale   up    to  ocean    waves  observed   from 
higher     altitudes. 

B.       LASER  RADAR  ALTIMETER  SYSTEM 

The  laser  system  used  was  a  working  model  of  a  laser 
radar  altimeter,  being  developed  by  a  companion  project. 
The  transmitter  consisted  of  a  0.905  micrometer  gallium 
arsenide  diode  laser  driven  by  a  pulser  having  a  pulse 
width  of  160  nanoseconds  at  half  height.  The  transmitter 
output  beam  divergence  was  1.4  degrees  by  2.8  degrees.  A 
silicon  avalanche  photodiode  detector  with  an  8  mm  focal 
length  lens  was  used  to  receive  the  reflected  laser  signal. 
The  aperture  field  of  view  was  circular  with  a  divergence 
angle  of  3.6  degrees,  totally  encompassing  the  area  illumi- 
nated by  the  transmitter.  A  10  nanometer  bandwidth  multi- 
layer film  filter  was  included  in  the  detector  optics  to 
eliminate  noise  due  to  sunlight,  for  daylight  operation. 
The  laser  system  was  mounted  on  an  arm  extending  3  feet 
from  the  bridge  railing.  There  were  no  bridge  supports  or 
obstructions  near  enough  to  produce  unwanted  reflections. 
The    laser   was    109    feet    above    the   water. 

The  detector  output  was  continuously  displayed  on  a 
portable  oscilloscope.  Because  the  pulse  length  of  the 
returned  laser  pulses  was  only  160  nanoseconds,  it  was 
necessary  to  stretch  the  received  pulses  electronically  and 
assemble  them  to  form  a  continuous-wave  analog  signal  to 
permit     recording.        By     the    Nyquist     theorem,      the    maximum 
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component  frequency  of  this  signal  is  limited  to  500  Hertz 
by  the  interpulse  interval  of  one  millisecond.  It  was  thus 
possible  to  record  this  signal  on  a  frequency-modulated 
analog  tape  recorder,  a  Hewlett-Packard  HP-3960A  recorder, 
which  had  a  signal  bandwidth  of  dc  to  20  kilohertz.  The 
signals  were  recorded  simultaneously  with  recording  of  the 
video  pictures  of  the  glitter  pattern.  Synchronization  of 
the  video  recordings  and  laser  return  signal  recordings  was 
accomplished  by  recording  the  same  voice  track  on  both 
recorders  . 

The  laser  return  signal  magnitude  was  later  data-pro- 
cessed in  the  laboratory  with  a  DATA  Precision  Corp.  DATA- 
6000  waveform  analyzer  to  yield  a  mean  laser  return  signal. 
As  the  laser  return  signals  were  played  back  from  the 
recorder  into  the  analyzer,  16  sequential  waveforms,  each 
consisting  of  512  data  points  taken  at  1  msec  intervals, 
were  averaged.  This  process  was  repeated  5  times  for  each 
run  and  these  were  averaged  to  obtain  a  mean  value.  Suc- 
cessive averagings  of  slightly  different  regions  within  the 
ten  seconds  over  which  the  video  signals  were  averaged 
yielded  a  standard  deviation  for  the  magnitude  of  the  mean 
s  i  gnal  . 

C.       LASER  REFERENCE   MEASUREMENTS 

An  effective  Lambertian  water  reflectance  signal  was 
obtained  by  means  of  measurements  of  the  laser  return 
signal     from    a    white    sheet    of    paper.        These    measurements 
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were  made  in  the  field  shortly  after  the  measurements  of 
water  reflectance  signal.  The  white  paper  sheet  was  lo- 
cated at  a  distance  of  39  feet.  At  this  distance  the 
transmitted  laser  pattern  was  still  fully  within  the  boun- 
daries of  the  sheet  of  paper.  The  same  paper  was  separate- 
ly calibrated  in  the  laboratory  and  found  to  reflect  0.73 
times  as  much  light  at  normal  incidence  as  a  perfect  Lam- 
bertian  surface.  This  value  is  representative  of  most 
white  paper  surfaces.  The  calibration  procedure  is  dis- 
cussed in  Appendix  B.  The  effective  Lambertian  water  re- 
turn signal  was  then  obtained  by  multiplying  the  signal 
from   the   white   sheet   of    paper  at    39   feet    by    the    ratios 

(.0204/0.73)  x  (39/109)2  . 
The  factor  (.0204/0.73)  takes  account  of  the  difference  in 
reflectance  of  the  white  paper  sheet  and  the  idealized 
Lambertian  water  surface.  The  factor  (39/109)2  takes  ac- 
count of  the  different  distances  of  the  white  sheet  and  the 
water,  39  and  109  feet  respectively.  The  inverse  square 
law  applies    to   these   distances. 

The  effective  reflectance  coefficient,  Peff»  as  defined 
in  section  II  E,  was  then  calculated  by  dividing  the  mean 
laser  water  return  signal  by  the  effective  Lambertian  water 
return    signal .    i .e 

Mean  laser  return  signal 


peff  " 


Idealized  Lambertian  water  return  signal 
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D.      GLITTER   PATTERN  MEASUREMENT 

Glitter  patterns  were  produced  by  illuminating  the 
water  surface  with  a  quartz-halogen  lamp  of  a  type  used  as 
a  light  source  in  overhead  classroom  projectors.  This  pro- 
duced a  nearly  uniformly  illuminated  patch  on  the  water 
with  a  total  beam  angular  divergence  of  about  60  degrees, 
considerably  wider  than  any  of  the  observed  glitter  pat- 
terns. The  uniformity  of  the  illumination  was  verified  in 
the  laboratory  in  advance.  The  lamp  was  rated  at  360  watts 
at  54  volts.  Power  for  this  lamp  was  provided  by  a  small 
rotary  converter  generating  nominal  120  volts,  60  Hertz 
power  up  to  500  watts,  and  driven  by  12  volts  dc  from  an 
automobile-type  storage  battery.  The  lamp  voltage  was  var- 
ied, as  needed,  by  means  of  a  variac.  The  video  camera 
was  an  RCA  type  TC2055C  with  a  Vidicon  tube.  It  was  oper- 
ated on  12  volts  dc  directly  from  a  second  automobile 
storage  battery.  The  camera  performance  was  very  stable  as 
it  generated  its  own  line-scan  frequency  with  a  quartz 
oscillator.  The  HP3960A  FM  data  recorder  for  the  laser 
radar  signals  also  operated  with  120  volt  60  Hertz  power, 
obtained  from  this  same  battery  through  a  frequency  con- 
trolled square  wave  chopper  converter.  The  video  images 
from  the  RCA  camera  were  recorded  on  a  portable  8mm  video 
tape  recorder  that  was  part  of  an  Olympus  VX-801-KU  8mm 
video  camera  system.  This  system  also  displayed  a  contin- 
uous    monitor     image    of     the    scene    as    viewed    by    the    RCA 
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camera.  This  allowed  the  brightness  level  of  the  illumina- 
tion to  be  adjusted  for  proper  recording  and  also  permitted 
alignment  of  the  camera  with  the  major  or  minor  axis  of  the 
elliptic  glitter  pattern.  The  8  mm  tape  system  was  very 
compact  and  operated  on  its  own  internal  batteries.  As 
measured  in  the  laboratory,  the  resolution  of  this  combined 
system  exceeded  that  of  the  usual  home  VHS  format  by  about 
20*. 

The  lamp  and  video  camera  were  mounted  on  an  arm  exten- 
ding from  the  bridge  railing  at  a  point  about  6  feet  from 
the  laser  altimeter.  The  video  camera  was  oriented  first 
along,  then  across,  the  wind  direction,  so  that  sequential- 
ly the  major  and  then  the  minor  axis  of  the  elliptical 
glitter  pattern  would  be  aligned  horizontally  in  the  re- 
corded video  image.  The  16  mm  camera  lens  was  focused  at 
infinity  with  an  aperture  setting  of  f/1.6  .  The  16  mm 
focal  length  gave  a  wide  enough  field  of  view  to  include 
the  complete  glitter  pattern  in  every  case.  Recording 
sequences  of  15  to  30  seconds  were  made  at  several  differ- 
ent illumination  levels  in  order  to  be  sure  to  avoid  satur- 
ating  the   camera. 

Later,  in  the  laboratory,  the  8  mm  video  tape  recor- 
dings were  copied  into  a  Sony  Superbeta  VCR  which  had  more 
flexible  reproducing  features  than  the  8  mm  field  unit. 
The  resolution  of  this  VCR  was  previously  measured  and 
found  to  exceed  that  of  the  usual  home  VHS  format  by  50*. 
No   detectable   degradation   occurred   in    this    copying    process. 


38 


The  output  of  the  video  recorder  was  sent  to  a  Tek- 
tronix 468  digital  oscilloscope.  Successive  single  hori- 
zontal TV  scan  lines,  lying  along,  or  perpendicular  to,  the 
glitter  pattern  major  axis,  were  displayed  on  the  scope 
using  the  "B  sweep  delayed"  mode.  The  waveform  displayed 
was  a  time-varying  signal  proportional  to  the  intensity 
distribution  along  the  major  or  minor  glitter  pattern  el- 
liptic axis.  The  oscilloscope  was  then  set  up  to  digitize 
and  average  256  such  sweeps  and  to  store  the  accumulated 
waveform.  The  real  time  period  represented  by  the  256 
sweeps  was  256/30  or  8.5  seconds.  The  computer  program  in 
appendix  A  was  written  to  allow  the  Tektronix  468  to  inter- 
face with  an  IBM  PC/XT  through  an  IEEE488  interface  bus 
[Ref.  7].  With  the  computer  as  controller  for  the  opera- 
tion, the  oscilloscope  digitized  and  transmitted  the  wave- 
form data  to  the  computer.  The  digital  information  was 
then   stored    in    memory   or   on   disk   for   further   processing. 

The  Tektronix  oscilloscope  was  then  adjusted  so  that  it 
displayed  a  TV  line  through  an  image  of  the  prevailing 
water  surface  outside  the  glitter  pattern.  Sweep  averaging 
and  transmission  of  a  waveform  to  the  computer  memory  were 
again  carried  out.  This  provided  a  background  signal  which 
was  subtracted  from  the  glitter  waveform  in  the  computer. 
The  computer  was  then  programmed  to  plot  the  glitter  pat- 
tern profile  on  an  Epson  FX80  dot  matrix  printer.  Figures 
3     and     4     show     the     resultant     intensity     profiles     along     the 
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major  and  minor  elliptic  axes  from  a  representative  run. 
Gram-Char  1  ier  standard  deviation  values  <ju  and  oQ  were  then 
measured  graphically  from  these  profiles  by  measuring  the 
half  width  at  a  height  equal  to  0.659/1.11  times  the  peak 
value,    as    discussed   earlier. 
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IV.    RESULTS 

The  measured  values  of  Peff»  the  ratio  of  observed 
laser  radar  return  to  that  expected  for  Lambertian  reflec- 
tion, were  plotted  as  a  function  of  the  corresponding 
values  of  1/(2<tu<tc),  in  Figure  5,  where  <ru  and  aQ  were  the 
measured  upwind  and  crosswind  Gram-Char  1 ier  standard  devia- 
tion values  for  the  glitter  patterns.  The  data  values  are 
tabulated  in  Appendix  C.  Use  of  the  function  l/(2<*ucrc)  for 
the  abscissa  was  suggested  by  the  model  discussed  in  sec- 
tion II.  This  functional  form  organized  the  experimental 
points  so  as  to  fall  fairly  well  along  a  straight  line. 
The  solid  straight  line  was  a  least  squares  solution  for 
the  data  points  taken  here,  shown  as  the  circles  in  that 
figure.  The  correlation  coefficient  was  r  =  0.96  for  the 
points  relative  to  this  line.  This  line  had  a  slope  of 
0.92   and    an    intercept    of    44    on    the   ordinate   axis. 

This  empirical  relationship  provided  a  basis  on  which 
to  predict  the  expected  reflectance  of  a  pencil  beam  from  a 
rough  water  surface  on  the  basis  of  only  the  remote  sensing 
measurement  of  the  glitter  pattern  profile.  This  made 
possible  the  field  evaluation  of  working  models  of  laser 
radar  altimeters  with  the  help  of  a  simple  technique.  The 
functional  relationship,  and  also  the  body  of  data  were  the 
first  of  this  sort  in  this  field.  Without  this  technique, 
the  field  performance  of  a  given  laser  altimeter  system 
would    require    expensive    concurrent    logistic    support    to 
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Figure  5.  Observed  effective  reflectance  relative 
to  Lambertian,  Peff>  as  a  function  of  1/(2  o  a Q)f 
plotted  as  circles.  The  solid  line  is  a  least 
squares  fit  of  a  straight  line  to  these  points.  Two 
data  points  plotted  as  crosses  are  based  on  Petri's 
data  [Ref.  2],  The  dashed  line 
model    proposed    in    this    thesis. 


represents    an    ideal 
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determine  the  roughness  state  of  the  water.  The  latter  was 
likely  to  be  inferred  from  meteorological  and  oceanographic 
data,  and  to  be  only  the  deduced  probable  steady  state 
condition,  whereas  the  technique  reported  here  directly 
measured    the    momentary    optical    state    of    the    ocean    surface. 

The  numerical  values  of  Peff  were  of  immediate  utility. 
The  lowest  value,  represented  by  the  intercept  on  the  y 
axis,  was  44.  Although  extrapolation  to  the  axis  had  con- 
siderable uncertainty,  it  did  imply  that  for  the  worst  case 
the  ratio  of  reflectance  to  that  for  Lambertian  was  con- 
siderably above  unity,  at  least  by  somewhat  over  one  order 
of  magnitude.  Without  evidence  such  as  this,  a  working 
altimeter  would  presumably  have  had  to  be  designed  to  cope 
with  a  case  where  a  unity  ratio  to  Lambertian  might  have 
occurred.  This  factor  of  44,  or  somewhat  over  one  order  of 
magnitude  increase  in  signal,  represented  a  large  reduction 
in  the  design  requirement  to  be  sure  of  a  return  signal. 
This  had  immediate  importance  in  the  ability  to  achieve  the 
design    requirements. 

That  the  observed  signals  were  directly  relatable  to 
reflection  from  a  Lambertian  surface  made  it  possible  to  do 
most  of  the  testing  work  in  the  laboratory.  Various  evolu- 
tions of  laser  radars  could  be  tested  with  a  Lambertian 
surface  such  as  a  white  sheet  of  paper.  The  signals  could 
then  be  scaled  to  those  for  a  water  surface  at  any  required 
distance.  The  range  of  field  signal  return  values  could 
then    be    evaluated    in    terms    of    Figure     5. 
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The  model    proposed   in    Section    1 1 B ,    which  predicted  that 

peff  =  t  Ho  ^  ^2<Tuffc^ »  reQuired  a  straight  line  through  the 
origin  with  slope  [  HQ  ]  =  1.11  in  Figure  5.  This  is  repre- 
sented by  the  dashed  line  in  that  Figure.  An  alternative 
prediction  is  for  the  slope  to  be.[HQ]/4  ,  depending  on  the 
view  taken  regarding  the  solid  angle  of  the  return  scatter- 
ing cone,  as  described  in  Appendix  D.  Although  neither 
alternative  represented  a  good  fit,  the  data  points  did  lie 
along  a  straight  line,  with  slope  not  far  from  1.11  .  The 
existence  of  a  theoretical  model  that  agreed  with  the 
observed  behavior  in  more  detail  would  have  been  desirable, 
but  it  was  not  necessary.  The  present  experiment  had 
achieved  the  stated  objective  of  providing  a  means  of 
determining  the  optical  state  of  the  rough  water  by  a 
remote  sensing  technique,  where  none  existed  before  this 
work.      The  result   was   semiempirical  ,    but   was   useful. 

The  fact  that  the  proposed  model  did  not  predict  the 
observed  behavior  could  be  viewed  as  an  interesting  chal- 
lenge. That  the  observed  effective  reflectance  was  larger 
than  expected,  by  somewhat  less  than  a  factor  of  four,  and 
that  the  straight  line  had  an  intercept,  indicated  that 
this  was  a  problem  unlikely  to  be  related  to  miscellaneous 
small  increases  or  losses  of  light.  Some  of  the  immediate, 
ad  hoc,  explanations  of  this  discrepancy  that  come  to  mind, 
seemed  to  predict  less  reflectance,  rather  than  more.  For 
example,  if  the  variation  of  specular  reflectance  with 
angle    were    invoked    for     light     reflected    at     large    angle    of 
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incidence  into  the  far  wings  of  the  Gram-Char  1 i er  distribu- 
tion, then  less  light  should  have  returned  along  the  axis, 
rather  than  more. 

Another  effect  producing  deviation  toward  too  little 
reflection  was  that  the  laser  "pencil"  beam  had  a  finite 
divergence  of  1.4  by  2.8  degrees.  Thus  the  reflectance  was 
really  the  average  over  the  curved  top  of  the  Gram-Char  1 i er 
distribution.  Thus  the  average  observed  would  have  been 
slightly  less  than  the  peak  value.  Correction  for  this  is 
planned  in  subsequent  work,  but  this  correction  will  pro- 
duce an  increase  in  the  peak  reflectance,  rather  than  a 
decrease,  as  needed  to  fit  the  model. 

A  possibility,  in  the  right  direction,  in  that  it  would 
have  tended  toward  producing  the  higher  than  expected  val- 
ues of  reflectance,  was  that  some  spray  might  have  begun  to 
develop  as  the  wind  increased.  This  would  have  correspond- 
ed to  the  region  of  small  values  of  1/(2ctu<jc).  The  drop- 
lets would  have  acted  as  ret roref lector s  and  increased  the 
returned  signal  along,  or  near,  the  laser  beam  axis.  This 
effect,  called  the  "glory",  is  commonly  seen  in  looking 
down  from  an  aircraft  in  the  direction  of  the  shadow  of  the 
aircraft.  However,  this  explanation  seemed  quite  unlikely 
in  view  of  the  low  sea  state  encountered  on  the  reservoir. 

As  to  the  major  factor  of  four,  it  appeared  that  there 
was  no  real  reason  to  omit  it.  The  belief  that  it  possibly 
lay  in  misinterpretation  of  the  solid  angle  of  the  return 
cone  of  light  had  its  only  real  motivation  in  that  the 
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alternative  choice  led  to  a  factor  of  four.  It  seemed  more 
likely  that  the  trouble  was  fundamental,  perhaps  lying  in 
the  assumptions  made  by  Cox  and  Munk,  and  others,  that  the 
light  was  reflected  by  infinitesimal  flat  facets,  rather 
than  by  a  statistical  average  over  focusing  curved  sur- 
faces, where  the  surfaces  fortuitously  curved  to  be  retro- 
reflectors  contributed  dis proport ional 1 y  high  return  sig- 
nal. Such  an  effect  seemed  likely  to  be  much  more  pro- 
nounced with  a  nearly  point  source  at  short  range,  as  used 
here,  than  for  the  case  of  the  sun  as  an  extended  source  at 
infinite  range,  as  used  by  Cox  and  Munk.  This  point  is 
expected  to  receive  more  attention  in  the  succeeding  work 
on    this    project. 
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V.    COMPARISON   WITH   OTHER   RESULTS 

Direct  comparison  of  the  results  reported  here  with 
those  of  others  was  possible  only  for  one  bit  of  frag- 
mentary data.  No  other  measurements  of  laser  reflectance  as 
a  function  of  glitter  pattern  width  were  known  to  the 
author,  except  for  a  single  diagram  in  the  article  by  Petri 
[Ref.  2],  The  values  of  glitter  pattern  width  were  not 
quoted  there,  as  the  article  was  primarily  concerned  with 
measuring  the  peak  reflectance  as  a  function  of  wind  velo- 
city. However,  a  few  sample  glitter  pattern  profiles  were 
shown  in  one  figure.  That  figure  is  reproduced  here  as 
Figure  6.  It  was  possible  to  measure  the  sigma  values 
graphically  from  these  curves,  although  the  curve  shapes 
were  poor.  Those  curves  were  obtained  by  a  laser  scanning 
in  a  vertical  plane,  with  the  plane  rotating  slowly  through 
360  degrees  about  a  vertical  axis.  The  scans  were  averaged 
so  that  no  information  was  retained  on  the  relative  magni- 
tude of  upwind  and  crosswind  sigma  values.  Two  of  these 
curves  yielded  sigma  values  that  could  be  related  to  the 
results  in  this  thesis.  The  others  were  too  poorly  re- 
solved or  they  corresponded  to  reflectivity  above  or  below 
the  range  measured  in  this  work.  Two  pairs  of  these  curves 
were  for  essentially  the  same  windspeed  so  that  an  estimate 
of  the  internal  consistency  could  be  made.  These  pairs 
differed  by  1296  and  1496  for  the  peak  value.  Points  for  the 
two    utilizable     widths     were     plotted     with    the    symbol     +    on 
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•Exceptions    (5  minute  averages 
in  one   azimuth  position) 

All   data  was  collected  with  a 
scanner  speed  of  16  scans/sec 
and  a   field  of  view  of  24 
mil  1 iradians 


27  Jun  -  3  mph   (estimated) 


28   Jun    -    7.4   mph 


29  Jun   -    9.6  mph 

4  Oct    -   9.6  mph 
29  Jun   -   10.4   mph 


-24.00  -16.00  -8.00  0.0  *8.00 

Degrees    from  the    vertical    (0°) 


♦16.00 


♦24.00 


Figure    6.      Normalized    average    return    signal    versus 
laser/receiver    angle    as    reported    by    Petri     [Ref.    2]. 
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Figure    5    of    this    thesis.     They    fell    close    to    the    data    points 
obtained    in    this    work. 

Although  the  results  obtained  in  the  work  reported  here 
were  intended  to  provide  a  means  of  evaluating  laser  radar 
altimeter  performance  without  the.  necessity  of  knowing  the 
wind  speed,  it  was  interesting  to  relate  these  results  to 
data  taken  by  others,  where  the  wind  speed  was  measured. 
This  provided  some  confirmation  of  the  validity  of  the 
results     reported    in    this     thesis. 

An  indirect  comparison,  where  an  inferred  wind  speed 
allowed  intercomparison  with  related  experimental  data  in 
the  literature,  was  possible  by  combining  the  results  of 
Petri  [Ref.  2]  with  the  results  of  Cox  and  Munk  [Ref.  3], 
Petri  reported  laser  reflectance  data  in  terms  of  wind 
speed  for  waves  under  the  Chesapeake  Bay  bridge  near  Anna- 
polis, Md.  Figure  7  reproduces  the  results  of  Petri  for 
laser  return  as  a  function  of  wind  speed,  with  his  results 
plotted  as  the  crosses.  No  theoretical  model  was  offered 
in     that     article. 

Cox  and  Munk  gave  a  relationship  between  their  measured 
wind  speed  and  the  mean  square  upwind  and  crosswind  sigma 
values  for  the  glitter  patterns,  where  W  was  the  wind  speed 
in   meters   per   second,   as 

oQ2    +    au2    =    .003    +    5.12    x    10"3    W    ±.004       . 
W    has    been   called    the    effective    wind    speed    here,    as    it    was 
only    a    means    to    cross-connect    to   the    data   of   Petri.       Using 
the     above     expression,     effective     wind    speeds     were     calculated 
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Figure  7.  Plot  of  observed  effective  reflectance  rela- 
tive to  Lambertian,  Peff>  versus  wind  speed.  Circles 
represent  reflectance  data  obtained  in  this  work  plot- 
ted versus  inferred  wind  speed.  Crosses  represent  ref- 
lectance data  collected  by  Petri  plotted  versus  mea- 
sured wind  speed  [Ref.  2]. 
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for  our  data.  The  corresponding  laser  returns  are  plotted 
as  the  circles  in  Figure  7.  It  can  be  seen  that  the  data 
points  from  both  sources  cluster  into  a  broad  band  in  the 
same     general     region. 

The  approximate  agreement  of  the  data  reported  here 
with  the  combined  results  of  Petri  and  Cox  and  Munk  gave 
some  credence  to  the  ability  to  relate  laser  radar  return 
to    the    width    of    glitter    patterns. 

It  should  be  pointed  out  again  that  the  purpose  of  this 
work  was  to  obtain  a  measure  of  expected  laser  radar  return 
signal  from  rough  water  by  remote  sensing  techniques  and 
without  a  knowledge  of  wind  speed.  The  effective  wind 
speed  deduced  from  the  expression  of  Cox  and  Munk  was 
useful  only  to  permit  the  cross-connection  to  Petri's  data. 
These  wind  speeds  should  be  considered  here  only  as  a  means 
of  obtaining  order  of  magnitude  confirmation  of  the  optical 
results.  Intercomparison  was  somewhat  doubtful  because  the 
Cox  and  Munk  data  were  for  open  ocean  waves,  with  wind 
speeds  measured  at  9  and  41  feet  above  the  water.  The 
Petri  results  were  for  waves  under  the  Chesapeake  Bay 
bridge  near  Annapolis,  Md,  with  wind  speed  measured  at  60 
feet  above  the  water.  In  spite  of  these  uncertainties  in 
wind  velocity  it  was  interesting  that  utilization  of  these 
results  yielded  approximate  agreement  between  the  data  re- 
ported   here    and    that    of    Petri. 
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VI.    CONCLUSION 

This  work  showed  the  feasibility  of  a  remote  sensing 
technique  for  determining  the  expected  magnitude  of  laser 
radar  return  from  a  rough  water  surface  by  measurement  of 
the  size  of  a  simultaneously  measured  elliptical  glitter 
pattern.  The  technique  did  not  require  any  additional 
knowledge  of  the  water  surface  beyond  optical  measurements 
made  at  the  location  of  the  laser  radar  altimeter.  No 
previous  such  technique  had  been  in  existence.  The  effec- 
tive laser  reflectance,  Peff>  proved  to  vary  approximately 
linearly  with  l/(2<*ucrc),  where  cru  and  o Q  are  the  Gram- 
Charlier  standard  deviations  for  the  upwind  and  crosswind 
glitter  pattern  intensity  profiles,  respectively.  Peff  is 
the  ratio  of  the  laser  radar  return  signal  to  the  laser 
radar  return  signal  to  be  expected  from  the  water  if  it 
were  to  act  as  a  Lambertian,  diffuse  reflecting,  surface 
but  with  an  overall  reflectance  coefficient  equal  to  that 
for  normal  incidence  for  water.  A  least  squares  fit  to  a 
straight  line  gave  a  slope  of  0.92  and  a  y  intercept  of  44, 
with   a   correlation    factor,    r    =    0.96    . 

Plotting  Peff  as  a  function  of  1/(2<*u<tc)  was  suggested 
by  the  model  proposed  in  this  thesis.  That  model  indicated 
that  a  straight  line  for  Peff  as  a  function  of  l/(2cru<*c) 
should  pass  through  the  origin  with  a  slope  of  [HQ]/4,  or 
perhaps  a  slope  of  [  HQ  ]  ,  where  [  HQ  ]  =  1.11  ,  under  a  dif- 
ferent,   controversial,    interpretation.      The   solution   to  the 
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discrepancy  between  the  observed  data  and  the  suggested 
model  is  a  problem  that  remains  for  further  work. 

That  the  straight  line  did  not  go  through  the  origin, 
but  instead  had  an  intercept  of  44,  provided  important 
information.  This  said  that  the  worst  case  of  very  rough 
water  should  provide  a  return  signal  44  times  the  magnitude 
of  that  from  a  water  Lambertian  surface.  Although  such 
extrapolation  completely  to  zero  is  risky,  it  implied  that 
for  very  rough  water  the  smallest  laser  radar  return  is 
somewhat  more  than  one  order  of  magnitude  larger  than  would 
have  had  to  be  assumed  in  the  absence  of  this  information. 
This  would  permit  the  design  of  a  laser  radar  altimeter 
with  an  order  of  magnitude  less  power  than  would  have  to  be 
assumed  otherwise. 

Finally,  the  existence  of  these  results  provides  a  firm 
relationship  between  the  reflectance  of  a  rough  sea  surface 
and  that  of  a  Lambertian  surface,  such  as  a  large  white 
surface  at  limited  range  in  the  laboratory.  This  reduces 
the  need  for  field  tests  of  various  working  models  of  laser 
radar  al t  imeters . 
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APPENDIX  A 


COMPUTER  PROGRAM 


10  DIM  Y(520) ,DTA(520) ,BACGND(520 
20  DEF  SEG  =  &HD000 

2  5  KEY  OFF 

3  0  INIT  =  0  'ASSIGN  OFFSET  ADDRESS  OF  INIT 

4  0  ADDR%  =21  'CONTROLLER  ADDRESS 

5  0  LEVEL%  =  0 

60  CALL  INIT(ADDR%,  1EVEL%) 

7  0  TRANSMIT=  3  : RECEIVE  =  6   'OFFSET  ADDRESSES 

8  0  CMD$     ="MLA  TALK  3"       'MY  LISTEN  ADDRESS 

9  0  'DEVICES'  TALK  ADDRESS 
100  CALL   TRANSMIT   (CMD$ , STATUS%) 

110  PRINT   "TRANSMIT  STATUS  =  ",STATUS% 

120  R$  =  SPACE$(255)  'SET  LENGTH  OF  R$>=  LENGTH 

130  Rl$  =  SPACE$(255):  R2$  =  SPACE$(255) 

14  0  PRINT  "Set  up  for  data  -  push  any  key  when  ready  " 

150  V$=INKEY$:iF  V$=""  THEN  150 

160  CLS  'OF  THE  INCOMING  DATA 

17  0  CALL  RECEIVE (R$,  LENGTH% , STATUS %) 

180  PRINT  "RECEIVE  STATUS  =  ",STATUS% 

19  0  CALL  RECEIVE  (Rl$ , LENGTH2% , STAT2%) 

200  PRINT  "RECEIVE  STAT2  =  ",STAT2% 

210  CALL  RECEIVE (R2$,LENGTH3%/STAT3%) 

220  PRINT  "RECEIVE  STAT3  ",STAT3% 

230  IF  STATUS%<>0  THEN  PRINT"RECEIVE  ERROR  TRY  AGAIN" ;GOTO 

100 
240  FOR  I  =  178  TO  255 
250  DTA(I-177)=ASC(MID$(r$,i) ) 
260  NEXT  I'270  PRINT  R$,R1$,R2$ 
280  FOR  1=1  TO  255 
290  DTA(I+78)  =  ASC (MID$ (Rl$ , iOO 

3  00  NEXT  I 

310  FOR  1=  1  TO  180 

320  DTA(333)=ASC(MID$(R2$,I) ) 

33  0  NEXT  I 

340  FOR  1=1  TO  520  STEP  10 

350  FOR  J  =  1  TO  10 

360  PRINT  USING  "     # # # " ; DTA (I+J-l) ; 

370  REM  LPRINT  USING"      ###" ;dta ( I+J-L) ;' 380  NEXT  J 

39  0  PRINT 

4  00  REM  LIST 
410  NEXT  I 

4  20  INPUT  "DO  YOU  WANT  THE  DATA  RECORDED?  Y/N" ,  ANS$ 

430  IF  ANS$  =  "n"  OR  ANS$  =  "N",  THEN  GOTO  580 

440  INPUT  "INPUT  FILE  NAME",  FILEN$ 

450  IPUT  "input  header" ,A$ 

4  60  OPEN  FILEN$  FOR  OUTPUT  AS  #1 
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47  0  PRINT  #1,A$ 

480  PRINT  #1,  MID$(R$,1,80) 

490  PRINT  #1,  MID$(R$, 81,80) 

500  PRINT  #1,  MID$(R$,161, 13) 

510  PRINT  #1,  DATE$,  TIME$ 

520  FOR  1=1  TO  520  STEP  10 

530  FOR  J=l  TO  10 

540  PRINT  #1,  USING  "     ###"; DTA (I+J-l)  ; 

550  NEXT  J 

560  PRINT  #1,"  " 

570  NEXT  I 

580  R$  =  SPACE$(255)  'SET  LENGTH  OF  R$>=LENGTH 

590  Rl$  =  SPACE$(255):  R2$  =  SPACE$(255) 

600  »OF  THE  INCOMING  DATA 

610  PRINT  "Set  up  for  background  -  push  any  key  when  ready  " 

620  V$=INKEY$:IF  V$=,ni  THEN  620 

63  0  CLS 

64  0  CALL  TRANSMIT (CMD$, STATUS %) 

650  PRINT  "TRANSMIT  STATUS  =  ", STATUS % 

660  CALL  RECEIVE (R$,  LENGTH% ,  STATUS%) 

670  PRINT  "RECEIVE  STATUS  =  ",  STATUS% 

680  CALL  RECEIVE (Rl$ , LENGTH2% , STAT2%) 

690  PRINT  "RECEIVE  STAT2  =  "/STAT2% 

700  CALL  RECEIVE (R2$ , LENGTH3% , STAT3%) 

710  PRINT  "RECEIVE  STAT3  =  ",STAT3% 

720  IF  STATUS%<>0  THEN  PRINT"ERROR  READING  DATA  " : GOTO  580 

722  IF  STAT2%<>0  THEN  PRINT " ERROR" : GOTO  580 

724  IF  STAT3%<>0  THEN  PRINT " ERROR" : GOTO  580 

730  FOR  I  =  178  TO  255 

74  0  BACGND(I-177)=ASC(MID$(R$,I)  ) 

750  NEXT  I 

760  PRINT  R$,R1$,R2$ 

770  FOR  1=1  TO  255 

780  BACGND(I+78)  =  ASC (MID$ (Rl$ , I ) ) 

790  NEXT  I 

800  FOR  I  =  1  TO  180 

810  BACGND(I+333)=ASC(MID$(R2$/I) ) 

82  0  NEXT  I 

830  FOR  1=1  TO  520  STEP  10 

840  FOR  J  =  1  TO  10 

850  PRINT  USING  "     ###"; BACGND (I+J-l) ; 

860  REM  LPRINT  USING  "     ###" ;bacgnd (I+J-l) ; 

870  NEXT  J 

88  0  PRINT 

890  REM  LIST 

900  NEXT  I 

910  INPUT  "DO  YOU  WANT  THE  DATA  RECORDED?  Y/N" ,  ANS$ 

920  IF  ANS$  =  '  n'  OR  ANS$  =  "N",  THEN  GOTO  1080 

930  INPUT  "INPUT  FILE  NAME",  FILEN$ 

940  INPUT  "input  header" ,A$ 

950  CLOSE  #1 

9  60  OPEN  FILEN$  FOR  OUTPUT  AS  # 1 
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970  PRINT  #1,A$ 

980  PRINT  #1,  MID$(R$,1,80) 

990  PRINT  #1,  MID$(R$,81,80) 

1000  PRINT  #1,  MID$(R$, 161,13) 

1010  PRINT  #1,  DATE$,TIME$ 

1020  FOR  1=1  TO  520  STEP  10 

1030  FOR  J=l  TO  10 

1040  PRINT  #1,  USING  "     ###"; BACGND (I+J-l) ; 

1050  NEXT  J 

1060  PRINT  #1,"  " 

1070  NEXT  I "1075  YMAX=-512 

1076  YMIN=512 

1080  FOR  1=1  TO  520 

1090  Y(I)=DTA(I)  -  BACGND(I) 

1091  IF  Y(I)<YMIN  THEN  YMIN=Y(I) 

1092  IF  Y(I)>YMAX  THEN  YMAX=Y(i) 
1100  NEXT  I 

1120  FOR  1=1  TO  520  STEP  10 

113  0  FOR  J=l  TO  10 

1140  PRINT  USING  "     ####"; Y ( I+J-l) ; 

1150  NEXT  J 

1160  PRINT 

1170  NEXT  I 

1180  XMIN=100 

1190  XMAX=364 

1500  YWMIN=YMIN-.2*(YMAX-YMIN) 

1510  YWMAX=YMAX+.15*(YMAX-YMIN) 

1520  XWMIN=XMIN-.2*(XMAX-XMIN) 

153  0  XWMAX=XMAX+.15*(XMAX-XMIN) 

1540  SCREEN  2,1 

1550  VIEW  (0,0)-(639,199) 

1560  WINDOW  (XWMIN,YWMIN)-(XWMAX,YWMAX) 

1580  YLAB$="RELATIVE  INTENSITY" 

159  0  XLAB$="POSITION  IN  GLITTER  PATTERN" 

1600  CLS 

1620  IXX=PMAP(XMAX,0) 

163  0  IXM=PMAP(XMIN,0) 

1640  IYX=PMAP(YMAX,1) 

1650  IYM=PMAP(YMIN,1) 

1670  DRAW"BM=IXM; , =IYM; " : DRAW"M=IXM; ,=IYX; :DRAW"M=IXX 

=IYX;" 
168  0  DRAW"M=IXX;,=IYM;":DRAW"M=IXM;,=IYM;" 

1681  LOCATE  23,  25 

1682  PRINT  XLAB$ 

17  00  FOR  X=XMIN  TO  XMAX 

1710  IX=PMAP(X,0) 

1720  IY=PMAP(y(X) ,1) 

17  3  0  IF  X=XMIN  THEN  DRAW"BM=IX; , =IY ; " 

17  4  0  DRAW"M=IX; ,=IY;" 

1750  NEXT  X 

1751  LOCATE  6 , 4 : PRINT  "RELATIVE" 
1790  WIDTH  "lptl:",255 
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1800  LPRINT  CHR$(27) "A"CHR$(7) 
2000  V$=INKEY$:IF  V$=IMI  THEN  2000 
2005  LPRINT  CHR$ (27) "A"CHR$ ( 12) 
2010  SCREEN  0 
3000  KEY  ON 
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APPENDIX    B 

LAMBERT I AN  REFLECTANCE  CALIBRATION 
Calibration  of  the  Lambertian  reflectance  coefficient 
for  the  white  paper  used  as  a  reference  in  the  field  was 
obtained  by  measuring  the  ratio  of  laser  radar  return 
signals  from  a  small  segment  of  area,  At  ,  of  the  white 
paper  target,  to  the  return  signals  for  specular  reflection 
from  a   plate  glass   window. 

The   Lambertian   reflectance,    r^,    of   the  the  white  paper 
is    given   by 


;  rs  PL  rl4 
_ 

4  At  Ps  «s2 


rT   = 


where 


rg    =   Reflectance  of   the   plate  glass   window,    for   two 
surfaces    =    .0826    for    n    =    1.51, 

P^  =   Laser   signal    from  a  small   area   Lambertian 
target , 

Ps    =   Laser    signal    from  a    plate   glass   window, 

R^   =  Distance    to    the    Lambertian    target, 

Rs    =   Distance    to    the    plate   glass    window,    and 

At   =  Area  of   the  white   Lambertian   target,    small 
enough    to    lie    within    the   central    flat 
illumination    region   of    the    laser    transmitted 
beam. 

The    calibration    described    above    was    independent    of    the 

Lambertian     properties    of    the    calibrating    target     of    white 

paper.        That    is,     it    did    not    matter    if    the    target    did    not 
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follow  the  Lambertian  behavior  of  brightness  being  propor- 
tional to  the  cosine  of  the  angle  from  the  normal.  This 
was  because  the  calibration  technique  involved  only  mea- 
surements made  with  the  laser  radar  and  hence  with  a  small 
acceptance  angle   near   normal   incidence. 

The  flatness  of  the  plate  glass  window  used  in  the 
calibration  was  measured  by  determining  the  position  of  the 
reflected  image  of  an  incandescent  lamp  source.  This  per- 
mitted calculating  the  radius  of  curvature.  The  glass  was 
then  treated  as  a  curved  mirror,  although  this  correction 
was  very  smal 1 . 

The  value  of  Lambertian  reflectance  obtained  was  0.73  . 
This  is  comparable  to  values  of  about  0.75  given  in  many 
tables  of  values  for  the  Lambertian  reflectance  of  white 
paper  . 
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APPENDIX  C 


TABLE  OF  DATA 

Glitter  Profile        Effective  Laser     Inferred 
Standard  Deviations        Reflectance     Wind  Speed 

(Peff) 

Predicted  Observed 


17.4  46.8  13.7 

20.5  47.6  7.18 
22.5  42.0  5.60 
71.4  90.5  2.41 
175  150  1.72 
142  182  2.90 
228  213  1.07 
194  217  0.92 

39.8  89.8  5.69 

52.9  92.7  4.04 
19.2  73.0  8.45 
34.8  88.3  4.80 
30.2  72.3  4.43 
37.0  86.9  3.58 
95.8  116  2.39 
12.7  43.5  10.4 
19.2  42.3  11.0 


Run 

# 

Upwi  nd 

Crosswi  nd 
(ac) 

1 

.340 

.094 

2 

.224 

.121 

3 

.175 

.141 

4 

.134 

.058 

5 

.122 

.026 

6 

.163 

.024 

7 

.081 

.030 

8 

.084 

.034 

9 

.218 

.064 

10 

.181 

.058 

11 

.241 

.120 

12 

.135 

.118 

13 

.167 

.110 

14 

.150 

.100 

15 

.126 

.046 

16 

.217 

.201 

17 

.301 

.096 

62 


APPENDIX  D 

LIGHT  SCATTERING  CONE 
In  calculating  the  expected  laser  reflectance  from  a 
rough  water  surface,  it  is  necessary  to  evaluate  the  size 
of  the  cone  into  which  the  surface  reflects  an  incident 
narrow  laser  beam.  However,  from  geometrical  optics,  this 
cone  does  not  have  the  same  included  solid  angle  as  the 
cone  measured  by  observing  the  glitter  pattern  obtained 
with  an  imager  located  near  a  point  source  illuminator.  In 
the  latter  case,  the  light  reflects  at  exactly  normal 
incidence  from  the  facets  on  the  water  surface,  and  the 
outer  limits  of  the  glitter  pattern  (or  sigma  points,  since 
a  definite  outer  limit  does  not  strictly  exist)  occur  at  a 
definite  angle  of  inclination  of  the  water  surface.  How- 
ever for  a  vertical  pencil  beam,  the  reflected  light  re- 
turns into  a  cone  that  has  twice  the  angular  spread  in  both 
x  and  y  directions,  and  hence  a  solid  angle  that  is  four 
times  that  of  the  glitter  pattern  observed  at  the  source 
with  a  point  source  illuminator.  This  is  because  the 
surface  which  the  pencil  beam  encounters  contains  facets 
that  are  inclined  at  all  the  angles  which  the  glitter 
pattern  indicates  are  in  existence,  if  averaged  over  a 
significant  period  of  time.  Thus,  for  example,  some  of  the 
pencil    beam    light    will    encounter    surface     facets    that    are 
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inclined  at  the  angle  corresponding  to  the  sigma  point  of 
the  glitter  cone.  Because  this  light  is  now  not  at  normal 
incidence  to  the  facet  it  is  encountering,  the  angle  of 
reflection  will  be  equal  to  the  angle  of  incidence  and  the 
total  angle  of  deviation  will  be  twice  the  angle  of  inci- 
dence. Since  the  angle,  at  the  glitter  pattern  position 
corresponding  to  the  sigma  point,  represents  the  angle  of 
inclination  of  a  facet  surface,  the  angle  of  deviation  from 
the  vertical  of  the  light  from  a  vertically  downward  pencil 
beam  will  be  at  an  angle  equal  to  twice  the  angle  for  the 
glitter  pattern  sigma  point.  Since  the  deflection  can  be 
in  either  the  x  or  y  direction,  the  solid  angle  of  the  cone 
into  which  the  light  from  the  pencil  beam  is  reflected 
would  be  expected  to  have  four  times  the  solid  angle  of  the 
glitter    pat  tern . 

Unfortunately,  the  results  obtained  with  the  above 
interpretation  lead  to  a  predicted  laser  return  that  is  too 
small  by  approximately  a  factor  of  four,  when  compared  to 
the  observed  laser  return  signal.  At  this  writing,  the 
explanation  of  this  anomaly  is  not  at  hand.  In  order  to 
lead  to  a  reasonable  fit  to  the  observed  laser  return,  it 
has  been  necessary  to  assume  that  the  cone  into  which  the 
reflected  light  radiates  is  equal  in  size  to  the  glitter 
pattern  cone. 
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